The Nature of Vibrations of a 2-D Disordered Lattice Model with Nano-scale Order 
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Motivated by the unusual soft mode anomalies in the relaxor ferroelectrics, we present a numerical 
study of vibrational excitations of a disordered lattice model with chemical order on the nano-scale. 
We find an Ioffe-Regel crossover in acoustic modes separating propagating low-frequency modes 
from the intermediate frequency diffusons, and localized modes at the high-frequency end of the 
spectrum. At a fixed degree of disorder, we find increasing uncertainty in the momentum of long 
wave-length optical modes with decrease in the size of nano-ordered regions. We determine the 
cause for this to be a strong mixing between acoustic and optical modes, believed to be central to 
the model of waterfall phonon anomalies. 
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Nano-scale structure of materials has interesting effects 
on their macroscopic behavior and become yet another 
control parameter in tuning material properties 0. For 
example, relaxor ferroelectrics, such as PbMgiNb^Oa 
(PMN)0 and PbSciNbi 3 (PSN)Q, have tremendous 
technological potential in applications based on solid 
state transducers^, due to their giant electromechani- 
cal and dielectric response. Randomness and nano-scale 
regions in the ordering of polar fields are known to be cen- 
tral to these interesting properties^ El- In contrast to 
conventional ferroelectrics such as PbTi03 with a sharp 
ferroelectric transition, relaxors exhibit a diffuse ferro- 
electric transition, with a broad peak in dielectric con- 
stant as a function of temperature. 

Ferroelectricity in conventional ferroelectric materials 
is essentially a phonon-related property characterized by 
a soft mode near the ferroelectric phase transition 0, EL 
The soft phonon anomalies in the context of relaxors |9j 
were observed as "waterfall phonons" recently using neu- 
tron scattering, with a speculation that these are linked 
with nano-polar regions. Further studies of the "water- 
fall" phenomena based on neutron scattering were ex- 
plained using a mode-coupling theory [Iollll| for the cou- 
pling between acoustic and optical modes. While the 
work of Hlinka et al argued that the wave-vector k 
of the waterfall phonon is unlikely to be related to the 
size of nano-polar regions, the origin of coupling between 
acoustic and optical modes is not known. 

A relaxor, for example PSN, PST etc., is a solid solu- 
tion with a degree of disorder arising from the chemical 
ordering of Sc and Nb (for PSN) and Sc and Ta (for 
PST) at the B sites of ABO3 perovskite lattice. Cor- 
relation between such disorder and relaxor properties of 
PSN and PST were shown experimentally |3l Il2j. A re- 
gion of the size of « 2-6 nm with a specific 
kind of chemical ordering defines a nano-ordered region 
(NOR) with intense local polar fields[6j, lira] . Through 
the coupling with ferroelectric dipoles, they give rise to 
nano-polar regions (NPR), speculated to be linked with 



phonon anomalies called "waterfall phonons" in relaxors. 
Precise knowledge of the nature of phonons in the pres- 
ence of chemically nano-ordered regions is essential to es- 
tablish their possible link with the "waterfall phonons" . 

In this letter, we determine the nature of vibra- 
tional excitations in a model system which consists of 
chemically-ordered nano-regions in a percolating disor- 
dered matrix. While the model used in our analysis 
is a great simplification of relaxors, it allows us to de- 
termine some general aspects of phonons in a nano- 
ordered system. In particular, we find localization of 
modes at high frequencies and an Ioffe-Regel crossover 
at low freq uenc ies, similar to that found in amorphous 
systems [l7l Ilil Q^j. Our findings about the nature of 
long- wavelength optical modes helps in uncovering possi- 
ble links between nano-scale order and waterfall phonons 
in terms of the mode coupling analysis of Hlinka et [l l| . 

For the ordered case, our model system consists of two 
types of atoms A and B occupying (0,0) a and (|, \)a 
sites of a square lattice respectively, a being the lattice 
constant. The interaction between two atoms displaced 
from their mean position by Ui and uj is given by 

Uij = fc L (b||.Ui)(b||.Uj) + fc T (b_ L .u 1 )(b_ L .u j ) (1) 

where b|| and bj^ are the unit vectors parallel (longitudi- 
nal) and perpendicular (transverse) to the line joining 
two atoms respectively, k^ and kx being the longitu- 
dinal and transverse force constants respectively. Only 
harmonic interactions up to the second nearest neigh- 
bors ((l,l)a) have been included, giving 10 distinct inter- 
atomic force constants k£p n , k££? nn , kf^ nn , /cf f 2nn , 
fcABrm an( j ^ WQ ons it e interactions kA and ks, where nn 
and 2nn stand for the nearest and next nearest neighbors 
respectively. Values of the force constants (Table I) have 
been chosen such that the acoustic sum rule is obeyed, 
there are stable modes throughout the Brillouin zone, 
and acoustic and optical branches are well-separated. 

Disorder is introduced by randomly selecting two dif- 
ferent types of atoms and interchanging their positions 
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making sure that a given atom is interchanged only once. 
A configuration with nano-scale order is obtained by ex- 
cluding such an interchange in the circular nano-ordered 
regions (NORs) with a fixed diameter of a few nanome- 
ters. Centres of these circular regions are taken ran- 
domly, but avoiding any overlap between them. The de- 
gree of disorder is quantified with a disorder parameter 
S defined as 



S 



No. of A-A or B-B nearest neighbors 
Total no. of nearest neighbor 



(2) 



For a fully ordered case, 5 = 0.0, whereas 5 is 0.5 for 
the completely disordered case. For a disordered con- 
figuration A and B occupy one of the two sites of the 
square lattice with equal probability, resulting in new 
types of neighboring pairs of atoms, such as (A-B)2nn and 
(A-B)3„„. The force constants for their interactions are 
taken as (k^ A ""+kf BnTl )/2.0 and (k£ 42n ™+kf B2,m )/2.0 
respectively, where X = L,T. Similarly, nearest neigh- 
bor atoms of the same type can possibly be separated 
by (0.5,0.5)a, for which the force constants are given by 
^disorder =2 ^™, where Y = A,B. The (Y-Y) ml 
and (Y-Y) 2rm interactions for the ordered case are now 
the (Y-Y) 2nn and (Y-Y) 3 „„ respectively. The acoustic 
sum rule is imposed by correcting the diagonal elements 
of the force constant matrix[2(J. Using this model, we 
study vibrational modes of configurations consisting of 
32 x 32 lattice sites (2048 atoms) with different number 
of nano-ordered regions and S ss 0.4. 



TABLE I: Force constants (in ar. units) for different interac- 
tions along different directions 

Force Constants Values Force Constants Values 
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The spectral density function (Sfc(w)) is calculated for 
different configurations, which allows us to relate to the 
constant energy scan and constant momentum transfer 
scan obtained in the neutron scattering experiments. We 
use the contour plot of Sfc(w) in visualizing the dispersion 
of phonon modes for systems with disorder. The spectral 
density function [2l| is given by 



Sk(w) 



■up) 



2N 



(3) 



where N is the total number of atoms. uP is the frequency 
of the j th eigenmode, \a J k \ is the Fourier transform of the 




FIG. 1: Phonon dispersion obtained from (a) the dynamical 
matrix of ordered system, and (b, c, d) contour plots of Sk(u) 
for the ordered system, a nano-ordered system with fourteen 
NORs each of radius ~ 2.2 a and a fully disordered system 
respectively. 



real space eigenvector of the j mode. Phonon disper- 
sion of the ordered system along the (10) direction ob- 
tained from dynamical matrix and from Sfc(w) are com- 
pared in in Fig. At X point, the longitudinal optical 
branch has the highest energy, whereas the transverse op- 
tical branch is degenerate with the longitudinal acoustic 
branch. Finite width of the peaks in the contour plot of 
Sk {ui) arise from the width of Gaussian used to approxi- 
mate the delta function in Eqn J2J. Based on a reasonable 
agreement between the two descriptions, we use Sk(u>)- 
based description in the study of phonons of systems with 
disorder. Our results for phonons of nano-ordered and 
disordered configurations (Fig. dc)-(d)) show that the 
introduction of disorder results in mixing of the phonon 
modes and subsequent broadening of peaks in Sfc(w) as 
a function of both k and u>, k is no longer a good quan- 
tum number to describe these modes. While new modes 
with intermediate energies appear near the T point, long 
wavelength acoustic modes for k < 0.22 are least affected. 
This is consistent with the error-bars on phonon measure- 
ments using inelastic neutron scattering technique 0- 

To probe extended versus localized nature of modes 
in the presence of disorder and nano-ordered regions we 
calculate inverse participation ratio (IPR): 



IPR(u)j) = N 



a: 



■J 14 



(S, 



(4) 



where a? are the eigenstates of the phonon modes. IPR 
is N (1) for a fully localized (extended) mode. For a con- 
figuration with fourteen nano-ordered regions, IPR as a 
function of energy (Fig. Ufa)) is found to be quite similar 
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FIG. 2: For the configuration with 14 NORs: (a) IPR as a 
function of energy and (b) acoustic mode frequency and its 
inverse decay time as a function of wavevector, point of their 
crossover is indicated by a circle. 



to that obtained for amorphous Si^^- The high energy 
modes, between 5.0 and 6.0 units of energy, are highly 
localized and are called locons, while the rest are called 
extendons|l8|. Some modes around e = 2.0 have inter- 
mediate values of IPR exhibiting moderate localization. 
We investigated these modes further by examining dis- 
tribution of number of atoms as a function of their con- 
tribution to the IPR and find that they are similar to the 
localized modes with the same IPR and 5 < e < 6. How- 
ever, the former have non-negligible character of long- 
wavelength modes (evident in Fig 1) and the latter do 
not. 

Extendons can be classified into two categories, the 
propagons (these propagate ballistically) and the diffu- 
sons. Wave vector k is a good quantum number for 
propagons, but not for diffusons. The Ioffc-Rcgcl (I-R) 
cross-over 01 separates propagons and diffusons as a 
function of frequency. In the frequency domain, the full 
width at half maxima (£) of a peak in the spectral density 
function is used to obtain the decay time r of a propa- 
gating mode due to scattering (t^ 1 — n( UJ (k)) 19\. The 
crossover point between oj vs. k and r _1 vs. k gives the 
location (k) of the I-R cross-over. For the configuration 
with fourteen nano-ordered regions (Fig. Hfb)), we find 
the I-R crossover at k = 0.22 r.l.u. and u> = 2.2. The 
vibrational modes with energy between 2.2 and 5.0 are 
thus diffusons. 

To further affirm this classification of modes, we stud- 
ied their level spacing distribution (LSD). LSD, with 
its foundation in the random matrix theory (RMT), is 
known to exhibit Wigner-Dyson statistics for diffusons 
and Poisson statistics for locons 0] . We calculated level 
spacings distribution using spectral unfolding, in which 
the original spectrum is mapped into an ordered set of 



numbers with uniform density of states. We find that 
the LSD for modes with 2.2 < e < 5.0 obeys Wigner- 
Dyson statistics whereas that for modes with e > 5.0 
obeys Poisson statistics with slight deviation near zero 
spacing. This supports our classification of these modes 
as diffusons and locons respectively. 

To explore mixing among phonon modes due to disor- 
der, we introduce an energy correlation function /(e, e ) 
as an overlap between the modes of the ordered system 
(|a„)) with those of the disordered system (|a^)): 

/(e,e') = Z m , n \(a°\a%)\ 2 5(e-e m )5(e - e n ). (5) 

It is based on the representation of modes of a disordered 
configuration as a linear combination of normal modes of 
a fully ordered system. The contour plots of /(e, e ) (Fig. 
EJl clearly show that the acoustic modes of a nano-ordered 
system in the IR-crossover region emerge from mixing of 
acoustic and optical modes of the ordered system. At the 
fixed degree of disorder, this mixing becomes stronger 
with the number of nano-ordered regions. Secondly, dif- 
fuson modes in the optical region of the spectrum of a 
nano-ordered system arise mainly from mixing among op- 
tical modes with different wave vector. We find that the 
effects of disorder become stronger as the number (size) 
of nano-ordered regions increases (decreases). 
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FIG. 3: Energy correlation between the modes of ordered 
system and that of a configuration having (a) 8 nano-ordered 
regions and (b) 14 nano-ordered regions. The straight line 
shows perfect correlation or NO mixing. 

To understand the nature of diffusons, we show in Fig. 
Hthe fc-dependence of S u (k) for uo ~ 4.292 ± 0.034 cor- 
responding to peaks of S at the Y point for a given con- 
figuration. It shows that the probability of finding a 
long wavelength optical mode decreases as the number 
of nano-ordered regions increases, as many of these de- 
grees of freedom precipitate into the acoustic branches 
due to strong mixing. Correspondingly, the uncertainty 
in momentum hk, which is given by the full width at half 
maximum (FWHM) of a peak in spectral-density also 
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In summary, we find that chemical disorder on a 2- 
D lattice results in phonon modes of type: propagons, 
diffusons and locons. At the fixed degree of disorder, 
reduction in the size of nano-ordered regions results in 
increase (decrease) in the width (height) of peaks cor- 
responding to long wavelength optical modes as a func- 
tion of k. This correlates with corresponding increase in 
the mixing between optical and acoustic modes. While 
the length-scale of chemical ordering needs to be large 
enough (a few nano-metre) for observing optical phonon 
modes with a definite momentum, it also has to be small 
to give adequate mixing with acoustic modes responsible 
for the "waterfall phonons" . 



FIG. 4: Probability of finding an optical mode in the long 
wavelength limit, "n" denotes the no. of nano-ordered regions 
present in the system and "r" is the radius of each nano- 
ordered region. These modes have energy oj=4.292±0.034 
units. 
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increases. For configurations having fourteen or more 
nano-ordered regions, there is no well-defined peak for 
which the FWHM can be estimated. 

Due to scattering from chemical disorder, optical 
modes (propagons) of the ordered system become diffu- 
sive. We suggest that the diffusons should be relevant 
to the diffuse dynamical properties of relaxors or dif- 
fuse phase transition behavior in ferroelectrics. While 
the phonons considered here are hard (high frequency) 
and treated within harmonic approximation, those in real 
ferroelectrics are soft and exhibit strong anharmonicity. 
Scattering of soft modes due to disorder is expected to be 
even stronger and their anharmonicity will lead to scat- 
tering among phonons. As the present work is on a 2-D 
system and some of its features are likely to be qualita- 
tively different in 3-D, precise implications for the real 
relaxors need further investigation. 

In connection with the waterfall phonon anomaly, our 
results show that nano-ordered regions of large enough 
size are essential for the observation of phonons with 
a definite wave-vector. Secondly, their size does af- 
fect the strength of mixing between optical and acous- 
tic modes, though not the wavevector corresponding to 
the I-R crossover. The former is the basis of mode- 
coupling theory and latter is consistent with the con- 
clusion of Hlinka et al0. We note that the NPR's (not 
NORs) are present in conventional ferroelectrics too (lik e 
the chain- like dynamic polar distortions in KNb03|22j) 
over a narrow temperature range above T c and "water- 
fall phenomenon" should be more general. As the acous- 
tic modes essentially span configurations with inhomoge- 
neous strain, the coupling of inhomogeneous strain with 
polar optical modes seems to be the common ori gin of 
"waterfall phonons" and also the relaxor dynamics 23] . 
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